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Section  1 
INTRODUCTION 


1 . 1  BACKGROUND 

Classical  Ekman  theory  considers  a  steady  wind 
stress  acting  over  a  homogeneous  ocean  of  infinite  depth 
and  constant  eddy  viscosity.  The  surface  current  deflects 
45°  to  the  right  of  the  wind  direction  (in  the  Northern 
Hemisphere),  and  the  deflection  angle  increases  with  depth. 
The  current  magnitude  is  greatest  at  the  surface  and 
decreases  exponentially  with  depth. 

If  the  wind  velocity  were  to  change  suddenly,  the 
currents  would  evolve  toward  their  new  equilibrium  profile. 
This  approach  to  the  new  equilibrium  would  be  oscillatory, 
with  inertial  frequency  f  *  2Qsind>,  where  ft  *  2n/day  is  the 
rotation  rate  and  4>  is  the  latitude.  These  inertial  oscil¬ 
lations  can  often  be  difficult  to  observe.  Inertial  oscil¬ 
lations  are  transient  in  nature,  and  can  require  long  time 
series  of  current  records  to  make  identification.  At 
latitudes  where  f  coincides  with  tidal  frequencies,  inertial 
oscillations  can  be  difficult  to  distinguish  from  tidal 
motions.  Examples  of  these  difficulties  are  shown  in 
Webster  (1968). 

There  are  two  components  of  the  near-inertial 
wave  field;  a  global  and  a  local  component.  The  global 
component  is  generated  by  a  variety  of  mechanisms  at  lower 
latitudes  as  super-inertial  internal  waves.  If  these  waves 
propagate  away  from  the  equator,  they  appear  as  inertial 
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waves  when  they  reach  their  turning  latitude.  As  they  can 
reflect  off  the  bottom  and  top  ocean  surfaces,  and  off  ir¬ 
regularities  in  the  Vaisala  frequency  profile,  they  may  be 
observed  to  propagate  upwards  as  well  as  downwards.  Since 
these  waves  are  observed  in  locations  that  are  far  from 
their  geographical  areas  of  generation,  they  probably  can 
only  be  described  in  a  statistical  sense.  Fu  (1980)  shows 
that  the  global  component  of  the  near-inertial  wave  field  in 
the  deep  ocean  (over  smooth  topography)  can  be  described 
successfully  with  a  universal  spectral  model. 

There  is  also  a  local  component  of  the  near- 
inertial  wave  field.  The  strongest  evidence  that  a  com¬ 
ponent  of  the  wave  field  is  locally  forced  is  the  obser¬ 
vation  that  the  dominant  direction  of  the  phase  propagation 
is  upward.  According  to  the  internal  wave  dispersion 
relationship,  this  corresponds  to  downward  energy  propa¬ 
gation.  The  implication  is  that  the  surface  layer  of  the 
ocean  acts  as  a  source  for  these  waves.  In  brief,  the 
hypothesized  mechanism  is  as  follows.  When  the  wind 
stress  field  is  inhomogeneous  in  space,  a  vertical  velocity 
field,  known  as  Ekman  pumping,  deforms  the  interface  between 
the  surface  layer  and  the  stratified  interior.  If  these 
vertical  motions  are  resonant  with  near-inertial  internal 
waves,  then  energy  can  propagate  away  from  the  interface, 
into  the  deep  interior. 

The  purpose  of  the  present  study  is  to  investigate 
the  locally  forced  component  of  the  near-inertial  internal 
wave  field.  In  particular,  we  are  interested  in  the  propa¬ 
gation  of  near-inertial  motions  and  vertical  shear  from  the 
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surface  layer  into  the  upper  oc^an  thermocline.  It  behooves 
us  to  separate  the  problem  conceptually  into  two  parts.  In 
the  first  part  we  examine  the  generation  of  inertial  oscil¬ 
lations  in  the  surface  layer.  Much  of  the  shear  across  the 
base  of  the  wind-mixed  layer  is  due  to  these  inertial 
oscillations,  and  this  shear  is  an  important  factor  in  the 
deepening  of  the  mixed  layer.  Inertial  shear  across  the 
mixed  layer/thermocline  interface  gives  rise  to  shear 
instabilities,  which  help  to  entrain  water  across  the 
interface  and  erode  the  thermocline  (Pollard  et  al.,  1973; 
Kruass,  1981).  Thus,  the  inertial  oscillations  in  the  mixed 
layer  are  important  in  their  own  right. 

The  second  part  of  our  problem  is  the  trans- 
ferral  of  near-inertial  energy  from  the  mixed  layer  into  the 
thermocline.  Near-inertial  kinetic  energy  is  not  as  strong 
here,  but  the  background  stratif ication  in  the  thermocline 
allows  the  energy  to  propagate  into  the  deeper  waters,  in 
the  form  of  free  internal  gravity  waves.  These  waves  are 
excited  through  resonances  between  the  wind  field,  Ekman 
pumping,  and  the  internal  wave  field. 

This  report  reviews  the  interactions  between  the 
wind  stress,  the  Ekman  layer,  and  internal  waves  in  the  in¬ 
terior.  The  purpose  is  not  to  relate  an  exhaustive  list  of 
the  models,  theories  and  observations  of  inertial  motions. 
Instead,  our  purpose  is  to  lay  a  basic  foundation,  with 
strong  emphasis  on  the  physical  dynamics  which  underly  iner¬ 
tial  motions.  In  Section  1.2  we  develop  the  general  set  of 
equations  which  describe  near-inertial  motions.  In  Section 
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2  we  look  into  the  time  dependent  response  of  the  wind- 
forced  Ekman  layer.  A  number  of  different  cases  and  con¬ 
ditions  are  examined.  We  show  how  the  attenuation  rate  of 
inertial  oscillations  may  be  related  (indirectly,  through 
the  eddy  viscosity)  to  the  background  stratification.  In 
Section  3  we  explore  two  classes  of  resonance  conditions. 
The  first  class  involves  the  resonance  between  the  wind 
field  and  the  Ekman  pumping.  This  sort  of  resonance  is 
characterized  by  an  oscillatory  wind  field  whose  frequency 
is  in  the  neighborhood  of  f.  The  second  class  is  the 
resonance  between  a  propagating  Ekman  pumping  field  and  the 
internal  wave  field.  This  resonance  activates  wave  packets 
whose  phase  velocity  is  equal  to  the  propagation  speed.  In 
Section  4  we  discuss  the  role  of  future  research  in  under¬ 
standing  inertial  oscillations  and  vertical  shear. 

Wherever  possible,  we  attempt  to  cast  light  on 
the  structure  of  vertical  shear.  In  our  data  analysis 
report  (Rubenstein  and  Newman,  1981),  we  showed  evidence 
that  vertical  shear  may  be  a  better  indicator  of  near- 
inertial  frequency  oscillations  than  the  current. 

1.2  SCALING  THE  EQUATIONS  OF  MOTION 

We  wish  to  develop  a  set  of  equations  that  can  be 
used  to  describe  various  aspects  of  near-inertial-frequency 
motions.  We  begin  with  the  momentum  equations  for  an 
incompressible  fluid  with  negligible  molecular  viscosity 

p Dtu  +  2p(fix  u)  =  -VP  -p£  -  F  ,  (1.1) 
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and  the  equations  of  mass  conservation  and  continuity 

Dt P  *  - Y  ,  V*u  =  0  .  Cl.2a,b) 

The  symbols  p  and  P  denote  density  and  pressure,  u  =  (u,v,w) 
is  the  velocity  vector,  and  £  is  gravitational  acceleration. 
F  and  y  are  the  divergences  of  turbulence  momentum  flux 
and  mass  flux,  respectively,  due  to  unresolved  motion.  The 
symbol  Dt  denotes  the  substantive  derivative  3t  +  u.V  , 
and  n  is  the  local  component  of  the  earth's  angular  rotation 
vector.  We  define  a  coordinate  system  with  z  positive 
upwards;  x  and  y  are  the  conventional  zonal  and  meridional 
directions.  The  mean  ocean  surface  lies  at  z  =  0,  and  the 
flat  bottom  at  z  =  -H. 

We  are  interested  in  the  small  fluctuating  motions 
that  perturb  a  stationary  reference  state.  This  reference 
state  is  defined  by  the  stratification  pr(z)  which  satisfies 
the  hydrostatic  relation 

dz?r  ~  —  gPr(^)  «  (1.3) 

The  total  density  and  pressure  are  given  by 

P(x,y,z,t)  =  Pr(z)  +  P’(x,y,z,t)  , 

p(x,y,z,t)  =  or(z)  +  o'(x,y,z,t)  ,  (1.4a,b) 

where  P'  and  p'  are  perturbations  from  the  reference  state. 
We  assume  that  these  density  perturbations  are  small; 

P  ’  «Pr  , 
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so  we  make  the  Boussinesq  approximation,  in  which  the 
effects  of  fluctuations  in  density  on  the  momentum  balance 
are  neglected,  except  in  the  buoyancy  term.  We  define 
buoyancy  to  be 


(1.5) 


where  ,o0  is  a  representative  value  for  the  density.  Then 
(1.1)  and  (1.2)  become 


Dtu  +  2(£x  u)  =  -  —V P'  +  bz  -  F/P0  (1.6) 

-  -  Po  — 

D-fb  +  wN2(z)  =  -B  ,  V*u  *  0  ,  (1.7a,b) 


where  B  is  the  divergence  of  buoyancy  flux,  z  is  the  unit 
vector  in  the  vertical  direction,  and  N(z)  is  the  buoyancy 
frequency  defined  by 


*2<‘>  *-©(£)  ;  <1-« 

These  equations  are  too  general  for  our  purposes, 
so  a  number  of  assumptions  allow  suitable  simplifications. 
We  are  interested  in  near-inertial  period  motions.  The  time 
scale  of  interest  is  T  =  v/Q.  We  assume  that  both  u 
and  p  vary  over  the  same  horizontal  length  scale,  L.  This 
length  scale  is  that  over  which  appreciable  variations  occur 
in  surface  wind  stress  or  surface  buoyancy  flux.  We  assume 
that  the  aspect  ratio  6  =  H/L  is  very  small;  g^OClO"2). 
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We  nondimensionalize  as  follows; 


(x,y) 

= 

(x*,y*)L 

(u,v) 

= 

(u*, v*)U 

z 

a 

z*H, 

w 

= 

w*W, 

t’ 

a 

t*f0“1, 

P' 

a 

P*P0. 

b 

a 

b*B0 , 

N 

a 

0 

ss 

* 

S5 

( 1 .9a-h) 


where  asterisks  indicate  nondimensional  quantities,  fQ  =  2Q, 
and  N0  is  a  representative  value  of  N(z).  We  substitute 
(1.9)  into  (1.7b)  to  obtain  (dropping  asterisks) 


I  <UX  +  V  +  I  “z  *  0 


(1.10) 


If  wz  is  to  be  the  same  order  of  magnitude  as  the  hori¬ 
zontal  velocity  divergence,  then 


W  =  6U  ; 


(1.11) 


the  ratio  of  vertical  to  horizontal  velocities  is  of  the 
order  of  the  aspect  ratio  or  smaller. 

Next  we  scale  the  buoyancy  equation,  (1.7a); 


’t  *  (551) 


B0(ubx  *  Tby  *  Wbz>  +  l— I 


*  * 
u  b 


[(u'b')  *  (v '  b' )  +  6“1(w'b')J  .  (1.12) 


1-7 


The  overbars  written  over  the  buoyancy  flux  represent 

averages  over  the  time  scale  f0_l.  The  primed  quantities 

are  the  nondimensional  buoyancy  and  velocity  components  that 

vary  appreciably  over  time  scales  small  compared  to  f0-l, 

*  * 

and  their  dimensional  magnitudes  are  b  and  u  .•  We  neglect 
the  horizontal  divergence  terms  of  the  buoyancy  flux  because 
they  are  of  order  6  smaller  than  the  vertical  divergence 
term.  We  are  interested  in  small  velocity  perturbations, 
such  that  the  Rossby  number  U/f0L  is  much  less  than 
unity.  We  therefore  neglect  the  nonlinear  advection 
terms.  We  require  that  in  the  absence  of  small  time  scale 
fluctuations,  b^  should  be  balanced  by  wN^ .  Therefore, 
B0  =  6UN02/fo>  and  (1.12)  becomes 


bt  * 


wN 


(*k) 


(1.13) 


We  apply  similar  arguments  to  the  horizontal 
momentum  equation  (1.6)  to  obtain 


ut  -  sin<}>v  =  -Px 


(u’w1  )z 


v  +  sin^u 


( 1 . 14a ,b) 


where  terms  of  order  <5  and  of  order  U/f0L  have  been  neg¬ 
lected,  and  PQ  =  p0Uf0L  has  been  asserted. 
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The  vertical  momentum  equation  deserves  special 
attention.  We  drop  the  nonlinear  advection  terms  and  the 
horizontal  flux  divergence  terms  to  get 


52wt  _  ($cos$v  =  -Pz 


(1.15) 


We  define  N0  =  f0/5  ;  then  the  balance  is  hydrostatic; 


0  =  -P  +  b 
z 


(1.16) 


The  vertical  momentum  flux  divergence  is  of  order 
6(u*2/HfQU)  with  respect  to  the  pressure  gradient  Pz,  while 
the  horizontal  flux  divergence  terms  in  (1.14)  are  of  order 
(u*2/Hf0U).  Thus,  the  vertical  flux  divergence  in  (1.15) 
is  neglected.  We  also  drop  the  w*  and  costf-v  terms,  because 
they  are  of  order  6 2  and  6.  Nearly  all  investigators  drop 
the  cos<}>v  term,  the  so  called  "traditional  approximation", 
but  sometimes  the  w-t  term  is  retained.  It  is  necessary 
to  retain  this  term  for  the  case  of  high  frequency  motions, 
where  T  <<  fQ_1.  But  (1.15)  shows  that  for  near- 
intertial-f requency  motions,  the  nonstationary  term  should 
be  dropped,  and  it  is  inconsistent  to  retain  w^  without 
also  retaining  the  vertical  Coriolis  force  term,  cos^v.  A 
number  of  investigators  make  this  mistake,  and  we  wish  to 
point  out  this  inconsistency. 
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Having  determined  which  terms  we  wish  to  retain, 
we  list  below  the  dimensional  equations  of  motion  for  easy 
reference; 


ut  -  f0sin$v  =  -px  -  (u'w')z 
vt  +  f0sin<j)u  «  -Py  -  ( v'  w'  )z 
0  =  **Pz  +  b 

bt  +  wN2(z)  =  -(w'b ' )z 

ux  +  Vy  +  wz  =  0  ,  (1.17a-e) 

where  p  =  P'/pq  has  been  substituted,  to  simplify  the 
appearance  of  (1.17). 
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Section  2 
THE  EKMAN  LAYER 


r 


In  this  section  we  review  the  transient  behavior 
of  inertial  motions  in  a  surface  layer.  We  begin  though, 
with  a  brief  discussion  of  the  steady  state  balance  between 
the  Coriolis  force  and  vertical  stress,  the  classical  Ekman 
spiral.  Then  we  discuss  the  sudden  onset  of  surface  stress 
over  an  initially  quiescent  ocean;  horizontal  currents 
accelerate,  oscillate  at  near- i ner t ial  frequencies,  and 
asymptotically  approach  the  Ekman  spiral  velocity  profile. 

In  the  way  we  formulate  momentum  diffusion,  the 
stratification  enters  our  problem  only  indirectly,  through 
the  eddy  viscosity  parameter.  We  note  in  Section  2.4,  that 
a  constant  eddy  viscosity  may  be  most  appropriate  in  a 
continuously  stratified  ocean,  where  the  scale  of  momentum- 
transferring  eddies  is  not  a  function  of  depth.  In  a 
homogeneous  surface  layer  the  stress  is  nearly  constant,  and 
the  eddy  viscosity  might  be  better  described  as  a  linear 
function  of  depth.  This  model  leads  to  better  agreement 
with  observations  of  attenuation  rates  of  inertial  oscil¬ 
lations.  We  conclude  this  section  with  a  discussion  of  the 
situation  where  an  unstratified  mixed  layer  terminates  in  a 
density  interface,  and  overlies  a  stable  diffusion-free 
layer.  We  briefly  describe  how  shear  across  the  interface 
controls  the  entrainment  rate  which,  in  turn,  by  its  control 
of  the  mixed  layer  depth,  influences  the  character  of  shear 
in  the  Ekman  layer. 


The  equations  of  motion  used  in  this  section  are 
developed  as  follows.  Consider  a  horizontally  homogeneous 
unstratified  ocean,  at  latitude  <j>  ,  with  constant  eddy 
viscosity  K.  Equations  (1.17)  reduce  to 

ut  -  fv  -  Kuzz  =  0 


vt  +  -  Kvzz  =  0  , 

(2 . la ,b) 

where  f  =  2ftsin  <P  is  the  Coriolis  parameter. 

conditions  are 

The  boundary 

Kuz  =  ip(t)  z  =  0  , 

u  0  z  -*■  -  ®  , 

(2 . lc, d ) 

where 

A 

Tp  is  an  imposed  surface  stress  vector. 

that 

We  define  a  complex  velocity  U 

(2.1)  can  be  combined  to  read 

=  u  +  i v ,  so 

Ut  +  ifU  -  KUZZ  =  0  , 

KUZ  =  T0(t)  z  =  0  , 

U  *+•  0  Z  -  oe  , 

(2 .2a-c) 

where 

t0  is  the  complex  surface  stress. 

2.1  THE  EKMAN  SPIRAL 

Before  considering  the  solution  to  this  problem, 
we  first  examine  the  steady  state  situation,  where  the 


Coriolis  force  is  balanced  by  friction  in  the  ocean  inte¬ 
rior.  This  problem  is  easily  solved  by  setting  U*  in 
(2.2)  equal  to  zero.  A  method  of  solution  is  derived 
in  Appendix  A.  We  rewrite  the  solution  (A. 10)  in  a  more 
convenient,  easily  recognizable  form; 

u(z)  =  V0  exp(-  tt|z  1 /D)cos( ^/4  -  tt (z  |  / D ) 

v(z)  =  VD  exp(-  Hz  I /D)sin(  ^/4  -  tt|z]/D).  (2.3a,b) 

Here  a  surface  stress  is  applied  in  the  y-direction,  V0 
is  the  magnitude  of  the  surface  current,  and  the  Ekman 
depth  is  given  by 

D  ■  Tr(2K/f)1/2  .  (2.4) 

Figure  2.1  shows  the  Ekman  spiral,  and  we  note  that  at  the 
surface,  the  direction  of  the  current  is  45°  to  the  right 
(clockwise)  of  the  wind.  If  we  integrate  (A. 10)  from 
z  =  -“to00,  then  the  resultant  vector  is 

OO 

Ro  =  f  U(z>dz  s  ~i^To/f  .  (2.5) 

—  CO 

This  result  states  that  the  total  mass  transport  in  a  water 
column  is  rotated  90°  clockwise  with  respect  to  the  direc¬ 
tion  of  the  surface  stress  vector. 

It  is  also  of  interest  to  note  that  the  shear, 
defined  by  S  *  (u22  +  v^)1/2,  can  be  computed  from  (2.3); 

/2ttV 

S  =  — p-2  exp(-TT  |zj /D)  .  (2.6) 
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Figure  2.1  Ekman  spiral  for  pure  wind-driven  currents 
projected  on  a  horizontal  plane.  V0  is 
the  surface  velocity  vector  deviated  by  45° 
to  the  right  of  the  wind  vector  (Northern 
Hemisphere).  [After  Ekman  (1905).] 


The  shear  is  simply  an  exponentially  decaying  function  of 
depth,  with  an  e-folding  length  scale  D/it. 

2.2  TIME  DEPENDENT  SOLUTION  -  INERTIAL  TRANSIENTS 


We  return  now  to  the  time  dependent  problem, 
described  by  equation  (2.2).  We  adapt  Gonella's  (1971) 
method  to  our  present  needs  by  considering  three  cases; 
I)  an  initial  velocity  UQ(z)  =  U(z,0)  and  vanishing  surface 
stress  T0(t)  =  0,  II)  a  surface  stress  t0(t)  =  t(0,t) 
and  vanishing  initial  velocity  U0(z)  *  0,  and  III)  surface 
stress  lasts  for  a  finite  duration.  We  outline  our  method 
of  solution  in  Appendix  A.  The  general  solution  is  the  sum 
of  cases  I  and  II; 


where 


U(z,t)  =  Uj ( z , t  )  *  U0(z)  +  U i ( z  ,  t ) *  2t0(t)  , 

(z)  (t) 


-if  t  9 

U]_  (z  ,t  )  =  Y(t)  — - - ,-75  exp  (-z^/4Kt)  ,  (2.7) 

(4nKt ^ 

and  Y(t)  is  the  Heaviside  step  function.  The  asterisks 
denote  convolution  operations,  the  first  with  respect  to 
z  and  the  second  with  respect  to  t. 


2.2.1  Case  I:  Initial  Velocity  Profile 

When  the  initial  complex  velocity  U0(z)  *  U(z,0) 
is  nonzero,  and  the  stress  x0(t)  is  set  equal  to  zero, 
then  the  first  part  of  (2.7)  is  the  solution  for  U(z,t). 
For  large  t,  the  solution  asymptotically  approaches 

-ift  2 

U(z,t)  - - — —  R  exp  ( -z  /  4Kt )  ,  (2.8) 

(4nKt)l/^  ° 
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where  R0  =  R(t  =  0)  is  the  initial  resultant  vector  defined 
by 


R(t ) 


j U(z,t) 


dz 


(2.9) 


—  00 


This  result,  as  pointed  out  by  Gonella,  shows  that  the 
velocity  asymptotically  becomes  independent  of  z,  as  t 
becomes  very  large.  The  amplitude  decays  as  t-1/2,  and  the 
limiting  period  is  2ir/f,  the  inertial  period.  The  current 
at  each  level  tends  to  come  into  phase  with  the  resultant  of 
the  initial  distribution,  R0.  At  a  given  depth  z,  depending 
on  whether  UD(z)  leads  or  lags  with  respect  to  R0 ,  the 
period  theoretically  tends  to  the  inertial  period  from  above 
or  from  below  (see  Figure  2.2).  Gonella  claims  that  this 
effect  might  to  some  extent  be  responsible  for  observations 
of  spectral  peaks  being  distributed  in  the  vicinity  of  the 
inertial  frequency;  f  +  2  to  4%.  According  to  (2.3),  an 
initial  angular  difference  of  tt/2  between  UQ(z)  and  R0 
becomes  zero  after  about  ten  inertial  periods,  and  gives  an 
average  period  of  2ir/f  +  2.5%. 

For  the  situation  where  a  constant  stress  suddenly 
shuts  off  at  t=0,  then  we  have,  in  effect,  a  case  where  the 
initial  velocity  UQ(z)  is  the  Ekman  spiral  (2.3).  We 
showed,  in  equation  (2.5),  that  the  initial  resultant  vector 
R0  is  rotated  90°  to  the  right  of  the  wind  direction. 
From  (2.3),  the  depth  z  at  which  UQ(z)  is  rotated  90°  to 
the  right  of  the  wind  is  -D/4.  For  a  typical  value  of  eddy 
viscosity  K  =  100  cm2/sec,  the  value  of  D  at  midlatitudes 
is  on  the  order  of  50  m.  If  the  steady  wind  stress  were  to 
terminate  suddenly,  then  the  asymptotic  evolution  of  the 
current  profile  for  large  t  approaches  that  given  by  (2.4). 


Figure  2 


y  Direction  of  northorn 


.2  Initial  velocity  distribution,  uQ(z).  If  the 
initial  velocity  is  either  leading  or  lagging 
the  resultant  velocity  R0,  it  will  tend  to 
come  into  phase  with  R0  through  the  action  of 
viscosity.  Note  that  the  period  is  less  (T  < 
Ti)  for  that  part  of  the  distribution  which 
lags  and  is  greater  (T  >  Ti)  for  that  part 
which  leads.  From  Gonella  (1971). 
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At  depths  shallower  than  D/4,  about  12.5  m,  the  velocity 
vector  should  begin  to  rotate  with  a  frequency  greater  than 
the  inertial  frequency,  while  at  depths  deeper  than  D/4,  the 
rotation  should  be  at  a  lower  frequency. 

Kundu  (1976)  computed  high  resolution  kinetic 
energy  spectra  from  current  meter  observations  taken  in 
shallow  water  (100  m  deep)  off  the  Oregon  coast.  In  order 
to  achieve  statistical  stability,  he  ensemble  averaged 
spectra  from  eleven  depths  ranging  from  3.4  to  95.9  m. 
Figure  2.3  shows  two  spectral  peaks,  one  above  and  one  below 
the  inertial  frequency.  Although  the  unstratified,  later¬ 
ally  unbounded,  infinitely  deep  ocean  model  does  not 
directly  apply  to  these  observations,  Figure  2.3  is  sug¬ 
gestive  of  the  theoretical  mechanism  being  described  here. 

2.2.2  Case  II;  Surface  Stress 

When  the  constant  surface  stress  tq  is  turned  on 
at  t  =  0  over  an  ocean  initially  at  rest,  the  solution  for 
U(z,t)  is  given  by  the  second  part  of  (2.7).  Performing 
the  convolution  with  respect  to  t,  the  solution  is 

r  -if0  o 

U(z,t)  -  2t  /  - - — exp(-z2/4Ke)  d6  .  (2.10) 

°JQ 

This  is  Fredholm's  solution  for  the  response  after  a  sudden 
onset  of  wind  stress.  Figure  2.4  shows  this  solution  as  a 
function  of  time,  at  the  surface  and  at  depths  z  =  -.5D,  -D, 
and  -2D,  where  D  is  the  Ekman  depth.  The  tick  marks  along 
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High-resolution  spectra  of  KE  (=  u2  +  v2)  , 
averaged  over  all  eleven  depths.  Bandwidth  = 
0.002  cph.  Note  the  shift  of  the  inertial  peak 
to  0.064  cph,  and  the  secondary  peak  at  0.058 
cph.  From  Kundu  (1976). 


/ 


/«0 
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Figure  2.4 


Hodographs  showing  the  development  of  a  pure 
drift  current  (infinitely  deep  ocean)  with  time 
at  different  depths.  The  depth,  z,  is  given  in 
terms  of  the  thickness  of  the  upper  layer  of 
frictional  influence,  D,  and  the  time  in 
pendulum  hours  after  the  start  of  a  constant 
wind.  [After  Ekman  (1905).] 


2-10 


I 

I 

the  spirals  indicate  periodic  time  intervals,  and  are 
labeled  in  units  of  pendulum  hours.  The  current  vectors 
rotate  around  the  steady  state  current  vector,  and  asymp¬ 
totically  approach  the  Ekman  spiral  solution  (2.3). 

2.2.3  Case  III:  Surface  Stress  for  a  Finite  Duration 


When  an  initially  quiescent  ocean  is  suddenly 
driven  by  a  surface  stress  that  lasts  for  a  duration  d, 
we  might  approach  the  problem  by  giving  the  wind  stress 
the  form 


t o ( t )  =  io[y(t)  -  Y(t-d)]  , 


(2.11) 


where  Y  is  the  Heaviside  step  function.  We  substitute 
(2.11)  into  the  second  part  of  (2.7)  to  get 


U( z , t )  =  2tq  J  exp( -z2 /4K0 )  c 
t-d 


-ife 


(  4ttK0  ) 


1/2 


d0 


(2.12) 


An  asymptotic  expansion  of  (2.12)  as  t  -*■  °°  may  be  written 


U(z,t) 


2t  e"iTr/2 

(1  -  e-ifd)  e"ift 


f(4TiKt) 


1/2 


(2.13) 


As  in  case  I,  the  period  approaches  the  inertial  period,  and 
the  amplitude  decays  as  t”1/2.  The  term  (1  -  e-ifd)  acts 
as  an  amplification  factor.  If  the  duration  d  is  an  integer 
multiple  n  of  the  inertial  period  (d  *  2irn/f),  then  this 
factor  vanishes,  and  the  oscillations  decay  rapidly  after 
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the  wind  terminates.  If  the  duration  is  a  half-integer 
multiple  of  the  inertial  period,  d  =  2  tt(  n  +  1  /  2) /f ,  then 
the  amplification  factor  is  maximized. 

An  instantaneous  change  in  the  wind  stress 
contains  Fourier  components  at  all  frequencies.  The  deeper 
water,  though,  is  shielded  from  stress  changes  on  small 
time  scales.  At  some  depth  z,  stress  changes  are  com¬ 
municated  on  a  time  scale  of  order  z2/k.  This  suggests 
that  in  a  time  dependent  Ekman  layer,  motions  are  in  the 
angular  frequency  range  f  +  Af,  where  Af  K/z^.  If  current 
measurements  show  energy  outside  this  frequency  range,  then 
these  motions  are  due  to  phenomena  other  than  diffusive 
penetration  of  surface  stress  changes. 

2.3  COMPARISON  WITH  OBSERVATIONS 

Gonella  (1971)  compared  selected  aspects  of  the 
time  dependent  solution,  described  in  Section  2.2,  with 
current  measurements  taken  in  the  Mediterranean  Sea.  Three 
periods  of  observations  corresponded  to  times  when  the 
atmospheric  conditions  were  favorable  for  comparisons  with 
theoretical  results.  During  each  of  these  periods,  a 
relatively  long  period  of  calm  weather  followed  storms  which 
could  be  well  modeled  by  an  impulsive  rectangle,  conforming 
to  Case  III. 

Preceding  two  of  these  measurement  periods,  the 
average  winds  were  of  comparable  intensity,  about  14  m/sec, 
but  lasted  different  durations.  Currents  measured  during 
July  1964  at  20  m  depth  were  about  twice  as  strong  as  those 
during  July  1968  at  10  m  depth.  The  first  measurements  were 
taken  after  a  storm  event  that  lasted  for  l£  inertial 
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periods,  while  the  latter  measurements  followed  a  storm 
event  that  lasted  2  inertial  periods.  As  discussed  in 
Section  2.2.3,  the  asymptotic  response  approached  long  after 
termination  of  the  wind  is  maximum  for  wind  events  lasting 
half-integer  multiples  of  the  inertial  period,  and  vanishes 
for  integer  multiple  durations.  The  observations  are 
consistent  with  this  model  result. 

Figure  2.5  shows  the  attenuation  of  inertial 
oscillations  during  the  three  observational  periods.  In 
each  case,  the  attenuation  rate  was  greater  than  the 
theoretical  t-1/2  rate.  Gonella  finds  the  observed  attenu¬ 
ation  rate  to  be  related  to  surface  layer  stratification. 
During  the  July  1964  and  1968  measurements,  the  strati¬ 
fication  was  continuous  down  to  30-35  m.  The  anomalous 
attenuation  rate  during  December  1964  was  measured  in  a 
surface  layer  that  was  homogeneous  down  to  about  70  m.  We 
refer  the  reader  to  Section  2.4  for  a  possible  explanation 
of  this  anomaly. 

2.4  DEPTH-DEPENDENT  EDDY  VISCOSITY 

Gonella  (1971)  made  a  cogent  comment  concerning 
Ekman  layer  models.  He  noted  that  although  the  Ekman's 
model  was  developed  for  a  homogeneous  ocean,  the  assumption 
of  a  uniform  eddy  viscosity  is  more  realistic  in  a  con¬ 
tinuously  stratified  ocean.  There  is  a  hint  in  Figure  2.5, 
that  Gonella's  model  may  be  better  applied  in  a  stratified 
rather  than  in  an  unstratified  surface  layer. 

Madsen  (1977)  addresses  this  issue  by  solving 
the  Ekman  layer  problem  with  a  depth  dependent  eddy  vis¬ 
cosity.  He  starts  with  a  generalization  of  equation  (2.2); 


Figure  2.5:  Attenuation  of  inertial  currents  as  a  function 
of  time.  The  time  origin  has  been  chosen  as 
the  beginning  of  the  wind  stress,  and  the  time 
is  expressed  in  units  of  inertial  period.  From 
Gonella  (1971).  During  July  1964  and  1968.  the 
stratification  was  continuous  down  to  30-35m. 
During  December  1964,  the  surface  layer  was 
homogeneous  down  to  70m. 


Ut  +  ifu 


( z )  U  z]  z  * 


(2.14) 


where  0  =  u  +  iv  is  the  complex  horizontal  velocity  vector. 
A  linear  form 

K(z)  .  ku*  |z|  (2.15) 

is  assumed  for  the  eddy  viscosity,  which  is  presumed  to  be 
appropriate  in  a  homogeneous  surface  layer.  Here  <  =  0.4 
is  von  Karman's  constant  and  u*  =  (  |t0|/p)1/2  is  a  friction 
velocity  based  on  a  representative  value  of  surface  stress 
magnitude  |  tD j .  The  boundary  conditions,  as  in  Section  2.1 
and  2.2,  are  that  stress  be  continuous  at  the  surface,  and 
that  U  ■*  0  as  z  ■*  -°°.  The  general  solution  for  an  initially 
quiescent  ocean  is  of  the  form 

t  (t)  ~ift 

U(z,t)  =  -2— •  *  exp(-!z|/KU#t)  ,  (2.16) 

*  (t) 

where  the  asterisk  denotes  a  convolution  with  respect  to  t. 
This  solution  is  of  a  form  analogous  to  (2.7)  for  the 
uniform  eddy  viscosity,  but  we  call  the  reader's  attention 
to  an  important  difference.  The  asymptotic  time  dependence 
here,  with  the  depth  dependent  eddy  viscosity,  is  propor¬ 
tional  to  t“l,  whereas  (2.7)  has  a  t-1/2  proportionality. 
Referring  to  Figure  2.5,  we  see  that  Madsen's  model  better 
predicts  the  attenuation  rate  for  the  December  1964  data 
set.  If  we  suppose  that  additional  processes  besides 
turbulent  mixing  are  simultaneously  contributing  to  the 
attenuation  of  inertial  currents,  then  we  conclude  that  the 
models  and  data  are  consistent:  In  a  neutrally  stratified 
surface  layer,  inertial  currents  decay  with  a  rate  at  least 
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as  fast  as  f1,  and  in  a  stably  stratified  surface  layer, 
the  decay  rate  is  at  least  as  fast  as  t-1/2.  An  example 
of  an  additional  process  which  could  contribute  to  atten¬ 
uation  is  internal  gravity  wave  dispersion,  discussed  in 
Section  4  of  this  report. 

Madsen  illustrates  that  his  model  is  different 
from  the  constant  eddy  viscosity  model  in  other  respects. 
Figure  2.6  compares  the  vertical  profile  of  velocity  in 
hodograph  form  with  the  classical  Ekman  spiral.  Two 
important  differences  appear.  First,  the  deflection  angle 
between  surface  wind  stress  and  surface  current  is  much 
smaller,  about  10°,  as  opposed  to  45°  for  the  Ekman  spiral. 
Madsen  points  out  that  his  result  is  consistent  with 
observations  of  oil  spill  trajectories,  both  real  and 
simulated,  which  show  deflection  angles  of  order  10°  or 
less.  Second,  there  is  a  more  rapid  decrease  and  rotation 
of  the  current  vector  with  depth.  The  vertical  shear 
of  the  steady  state  current,  in  the  limit  of  small  depth 
(  z  <<  u*/f)  is 

u* 

S  «  (uz2  +  vz2)l/2sr^j  .  (2.17) 

Unlike  the  exponential  decay  of  the  Ekman  spiral  in  (2.6), 
the  shear  in  the  depth  dependent  eddy  viscosity  model  is 
approximately  inversely  proportional  to  depth.  The  shear  at 
z  =  0  in  (2.17)  is  infinite,  a  result  of  the  fact  that  since 
K(z)  vanishes  at  z  »  0,  (2.16)  is  a  divergent  integral  at 
the  surface.  Madsen  copes  with  this  problem  by  translating 
the  vertical  coordinate  by  a  distance  ks/30,  where  ks  is  a 
typical  value  of  surface  roughness,  about  4  cm  for  moderate 
(3-10  m/sec)  surface  wind  speeds. 
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Figure  2. 


:  Vertical  velocity  structure  of  a  pure  drift 

current  in  an  infinitely  deep  homogeneous  ocean 
of  infinite  lateral  extent,  comparing  the 
turbulent  Ekman  spiral  (•)  and  the  classical 
Ekman  spiral  (+).  From  Madsen  (1977). 
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The  temporal  evolution  of  the  depth-dependent 
eddy  viscosity  model  is  much  more  rapid  than  Fredholm’s 
solution  (2.10).  Figure  2.7  compares  the  buildup  of  the 
surface  current  after  the  sudden  onset  of  surface  stress. 
The  time  increments  are  labeled  in  pendulum  hours.  A 
pendulum  hour  is  equal  to  a  sidereal  hour  divided  by  the 
sine  of  the  latitude.  A  requirement  that  surface  current 
magnitudes  be  equal  was  imposed  on  the  solutions.  The 
solution  to  Madsen's  model  approaches  a  steady  state  within 
about  three  pendulum  hours.  The  reason  for  the  faster 
adjustment  rate  is  that  the  equivalent  value  of  eddy 
viscosity,  evaluated  at  z  =  -u*/2f,  is  about  50  times 
the  value  in  the  classical,  constant  eddy  viscosity  model. 

2.5  DENSITY  INTERFACE  LIMITED  EKMAN  LAYER 

So  far  we  have  discussed  the  behavior  of  an  Ekman 
layer  in  a  homogeneous,  infinitely  deep  ocean.  With  a 
sudden  change  in  surface  stress,  a  momentum  impulse  diffuses 
downward.  The  leading  edge  of  such  an  impulse  lies  at  a 
depth  of  order  (2Kt)l/2.  After  a  while,  the  Coriolis  force 
leads  to  a  cross  flow  which  soon  balances  the  stress.  As 
a  result  of  this  balance,  most  of  the  vertical  shear  is 
confined  within  the  Ekman  depth,  D  =  ir(2K/f  )  1 /2  . 

When  a  mixed  layer  of  depth  h  lies  above  a  stably 
stratified  layer,  where  turbulent  diffusion  is  suppressed  by 
the  buoyancy  force,  the  vertical  growth  of  the  surface  shear 
layer  is  stopped.  Csanady  and  Shaw  (1980)  developed  a  model 
that  describes  the  interaction  between  these  two  depth 
limiting  mechanisms.  They  suppose  that  entrainment  of  fluid 
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Figure  2. 


The  temporal  development  of  the  surface  drift 
current  due  to  a  suddenly  applied  uniform 
surface  shear  stress.  Time  from  time  of 
application  is  indicated  in  pendulum  hours. 
The  present  eddy  viscosity  model  is  shown  by 
solid  circles  and  Fredholm’s  classical  solution 
by  plus  signs.  From  Madsen  (1977). 
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from  the  stable  layer,  across  the  density  gradient  inter¬ 
face,  into  the  mixed  layer  can  be  related  to  the  stress  at 
the  interface.  If  the  velocity  in  the  stable  layer  just 
below  the  interface  is  zero,  then  the  interface  stress  may 
be  written  in  terms  of  the  eddy  viscosity: 

Kuz  =  CduCu^+v^) 1/2  at  z  =  -h  .  (2.18) 

where  c<j  is  a  drag  coefficient. 

The  entrainment  velocity  we  is  presumably  proportional  to 
some  power  of  the  interfacial  friction  velocity 

u**  =  [cd(u2+v2)J1/2  ,  (2.19) 

and  is  consequently  proportional  to  the  magnitude  of  the 
shear  Au  across  the  interface.  Csanady  and  Shaw  choose  a 
power  law  found  by  Csanady  (1978); 


w 

e 


(2.20) 


where  t  is  kinematic  viscosity  and  Ap  is  the  density  jump 
across  the  interface. 

After  an  initial  momentum  impulse  reaches  the 
mixed  layer  depth,  entrainment  begins.  As  the  mixed  layer 
depth  increases,  the  interface  shear  Au  decreases,  and  the 
entrainment  rate  also  decreases.  If  the  mixed  layer  depth 
is  much  greater  than  the  Ekman  depth  (h>>D),  then  negligible 
velocities  at  z  «  -h  give  rise  to  no  further  entrainment. 
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It  seems  likely  that  an  initially  shallow  mixed  layer  depth 
should  asymptotically  approach  a  small  multiple  of  the  Ekman 
depth. 


Csanady  and  Shaw  (1980)  numerically  integrated 
equations  (2.1)  but  rather  than  using  the  boundary  condition 
(2.  Id)  that  velocity  vanish  at  infinity,  they  imposed  the 
boundary  condition  (2.18).  The  surface  stress  was  taken  to 
be  in  the  y-direction,  and  the  initial  value  for  h  was  9  m. 
A  uniform  eddy  viscosity  was  parameterized  as 

K  =  u*h/R  ,  (2.21) 

where  u*  was  the  surface  friction  velocity  and  R  =  12  was 
a  Reynolds  number.  Sample  solutions  are  shown  in  Figures 
2.8  and  2.9.  At  any  given  depth,  the  velocity  is  a  super¬ 
position  of  a  constant  component  plus  an  inertially  oscil¬ 
lating  component.  When  the  components  at  the  mixed  layer 
bottom  are  in  phase,  the  entrainment  rate  quickens.  The 
inertially  oscillating  component  damps  more  quickly  in 
Figure  2.9  because  the  initial  mixed  layer  depth  is  shal¬ 
lower  . 

Csanady  and  Shaw  also  note  that  in  Figure  2.9. 
where  the  ratio  h/D  is  quite  small,  the  velocity  component  u 
(perpendicular  to  the  wind)  is  nearly  uniform  in  the  mixed 
layer,  approximately  u*2/fh.  The  implication  is  that  the 
interior  stress  varies  almost  linearly  from  a  maximum  at  the 
surface,  to  zero  at  the  interface.  The  wind  stress  is  dis¬ 
tributed  almost  evenly  across  the  mixed  layer  and  balanced 
by  the  Coriolis  force  of  a  nearly  constant  drift  to  the 
right  of  the  wind. 
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h  in  meters  VELOCITY  (cm/sec) 


u  (cm/stc) 


MIOOEPTH 


v  (cm/stc/ 


BOTTOM 


u  (cm/stc) 


Hodographs  at  (a;  surface,  (b)  middepth,  and  (c)  bottom 
mixed  layer,  with  R  =  12,  c^  =  3  x  10-4,  v*  =  1,  and  init 
depth  6m.  Fast  convergence  to  steady  state  values  is  par 
due  to  the  artificial  damping  associated  with  the  numeri 
scheme.  From  Csanady  and  Shaw  (1980). 
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Section  3 

NEAR-INERTIAL  INTERNAL  WAVES  IN  AN  F-PLANE 

When  a  surface  wind  field  is  spatially  inhomo¬ 
geneous,  it  generates  a  vertical  flow  velocity  field,  known 
as  Ekman  pumping,  at  the  base  of  the  surface  layer.  If  the 
wind  field  varies  in  time,  then  under  certain  conditions  the 
Ekman  pumping  couples  the  surface  layer  with  the  inviscid 
interior.  This  coupling  is  achieved  through  resonances  with 
free  internal  waves.  These  free  internal  waves  can  then 
propagate  into  the  deep  interior. 

In  this  section  we  discuss  a  hierarchy  of  reso¬ 
nances  and  models.  In  Section  3.1  we  discuss  the  resonances 
that  occur  between  the  wind  field  and  the  Ekman  pumping.  We 
show  that  a  rotational,  divergent  wind  field  that  oscillates 
with  frequency  has  two  resonances,  at  frequencies 
(f2+c2)1/2,  wh  ere  c  is  a  damping  rate  coefficient.  This 
damping  results  from  stresses  at  the  base  of  the  surface 
layer . 

In  Section  3.2  we  discuss  the  condition  for 
resonance  between  the  Ekman  pumping  and  internal  waves.  We 
show  that  a  wave  whose  phase  speed  is  equal  to  the  migration 
speed  of  the  Ekman  pumping  velocity  field  satisfies  the 
resonance  condition. 

We  are  not  interested  in  the  Ekman  pumping 
velocity  field  per  se,  but  in  the  vertical  shear  associated 
with  near-inertial  internal  waves.  In  order  to  describe  the 
locally  generated  internal  wave  field  in  terms  of  a  given 
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wind  field,  we  combine  the  concepts  developed  in  Sections 
3.1  and  3.2.  In  Section  3.3  we  discuss  Krauss'  (1976a, b; 
1978a, b)  model  of  resonance  directly  between  the  surface 
stress  and  the  internal  wave  field.  We  show  that  if  a  wind 
field  contains  energy  in  certain  horizontal  wavenumber 
ranges,  free  internal  waves  are  excited. 

Finally  in  Section  3.4  we  discuss  the  propagation 
of  internal  waves  and  their  associated  vertical  structure. 

3.1  EKMAN  PUMPING  IN  A  VISCOUS  SURFACE  LAYER 

We  begin  with  equations  (1.18),  simplified  here 
for  the  case  of  a  homogeneous  fluid; 

ut  -  fv  =  -3Z  u'w' 

vt  +  fu  =  -3Z  v'w' 

ux  +  Vy  +  wz  =  0  (3.1a-c) 

We  are  concerned  primarily  with  the  vertical  velocity  w^  at 
the  bottom  of  the  surface  layer.  We  integrate  (3.1)  over 
the  depth  h,  first  by  defining 

0 

(u,v,wh,T0x, t0Y)  =  |dz(u,v , 9zw, 9zu’ w’ , 9zv( w’ ) 

-h  (3’2) 

Here,  u  and  v  have  dimensions  1 e n g t h 2 / t i me ,  and 
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_t0  =  (T0x,TQy)  is  the  surface  stress  vector.  Equations 
(3.1)  become 

ut  -  fv  =  (u'w' )h  -  T0X 

vt  +  fu  =  (v'w')h  -  t  Qy 

ux  +  Vy  »  -wjj  .  (3.3a-c) 

We  propose  a  simple  bottom  drag  parameterization  for  the 
stresses  at  depth  h; 

(u’w’  )jj  =  -cu 

(TV)h  =  -cv  .  (3.4a,b) 

As  discussed  in  Section  2.5,  if  the  depth  h  is  chosen  to  be 
somewhat  greater  than  the  depth  u*/f  =  td then  these 
stresses  in  (3.4)  are  very  nearly  zero.  We  can  combine 
(3.3a, b)  and  (3.4)  to  yield 

v  +  2cv  +  (c2+f2)v  =  fT0*  -  cx0y  -  r0y  .  (3.5) 

Here,  the  dots  indicate  partial  time  derivatives.  The 
importance  of  this  equation  lies  in  the  fact  that,  if 
io  is  horizontally  homogeneous,  then  so  too  is  v.  The 
equation  in  terms  of  u  is  analogous.  If  both  u  and  v  are 
horizontally  homogeneous,  then  no  Ekman  pumping  occurs; 
wh  ■  0  by  (3.3c) . 

We  could  solve  (3.5)  and  examine  its  solution, 
but  we  are  ultimately  interested  in  the  internal  waves 


in  the  deeper  ocean  which  are  driven  by  w^.  Therefore  we 
take  the  derivative  of  (3.5)  with  respect  to  y,  the  deri¬ 
vative  of  the  analogous  u  equation  with  respect  to  x,  and 
with  the  help  of  (3.3c)  we  obtain 

wjj  +  2cwh  +  (c2+f2)wjj  = 

fz  • (yXTp)  +  C7*To  +  7‘— o  -  (3.6) 

where  z  is  the  unit  vector  in  the  vertical  direction.  In 
the  steady  state  case  with  no  damping,  (3.6)  simplifies  to 
the  well  known  result  for  Ekman  pumping; 

wh  =  z*Vx£0/f  . 

Because  of  their  similarity,  some  of  the  following  remarks 
concerning  (3.6)  apply  also  to  (3.5). 

Equation  (3.6)  describes  the  motion  of  a  damped 
harmonic  oscillator.  An  equation  of  this  type  has  two 
solutions;  a  complementary  solution  to  the  homogeneous 
equation  which  describes  the  free,  undriven  response,  and  a 
particular  solution  to  the  inhomogeneous  equation  which 
describes  the  driven  response. 

The  complementary  solution  has  the  general  form 

whc  *  (Acosft  +  Bsinft)e-Ct  .  (3.7) 

The  natural  frequency  of  the  free  oscillations  is  f,  and  the 
solution  is  damped  by  an  exponential  envelope.  Damped 
oscillators  usually  have  a  natural  frequency  which  is 
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dependent  on  the  damping  constant,  but  for  the  special 
case  of  (3.5)  and  (3.6),  this  dependence  exactly  cancels! 

The  particular  solution  has  the  form 

whp  .  f-1  sinft  e-ct  *  [7. t0  +  cV*lo] 

+  sinft  e~ct  *  pa-vxlo^  ,  (3.8) 

where  asterisks  denote  convolution  operations  with  respect 
to  time.  This  solution  is  obtained  using  the  method 
of  variation  of  parameters,  along  with  trigonometric 
identities.  In  the  case  of  no  damping  c  is  zero,  and  (3.8) 
can  be  simplified  using  the  method  of  integration  by  parts; 

Whp  *  sinft  *  z  •  (Vxx0)  -  cosft  *  V-Jo  .  (3.9) 

We  summarize  some  important  results  concerning  (3.8)  and 
(3.9) . 

1)  The  convolution  operations  imply  that  the  system  has  a 
memory  of  the  historical  surface  stress.  In  (3.8),  the 
e-ct  factor  produces  a  fading  memory. 

2)  The  divergence  and  curl  terms  indicate  that  the  Ekman 
suction  velocity  wjj  is  produced  by  horizontal  gradients 
in  the  surface  stress. 

3)  The  Ekman  suction  velocity  is  modulated  by  inertial 
oscillations,  but  is  not  a  function  of  the  depth  h. 
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We  can  look  at  the  particular  solution  for  two 
selectively  chosen  surface  stresses. 


3.1.1  Rotational,  Solenoidal  Wind  Field 

We  consider  only  the  curl  component  of  the  wind 
stress.  When  jrQ  oscillates  with  a  frequency  u.'0,  we  might 
express  it  in  the  form 


_Io  ■  CO,  <J>0 ( x )  cosuot)  . 

The  particular  solution  to  (3.6)  may  be  written 

f  "  <$ ) 

%P  =  -  f2  -  c2)2“  +  4.V]1/2 


(3.10) 


6 


(3 . 1  la , b) 


The  damping  gives  rise  to  a  phase  shifts.  The  resonant 
frequency  is  found  by  differentiating  the  amplitude  with 
respect  to  u)Q,  and  equating  to  zero; 


wr  =  (f2  _  c2)1^2  . 


(3.12) 


The  effect  of  damping  is  to  reduce  the  resonant  frequency 
below  the  inertial  frequency. 


It  is  also  interesting  to  examine  the  complete 
solution  for  wjj  in  the  case  of  no  damping  (c  •  0). 
If  (3.10)  describes  a  surface  stress  that  turns  on  sud¬ 
denly  at  t  *  0,  then  appropriate  initial  conditions  are 
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wh(0)  =  wh(0)  =  0.  The  complete  solution -is  a  combination 
of  complementary  a..d  particular  solutions,  and  can  be 
expressed 


2f  Mo  ( w  +  f)t  (ui_  -  f)t 
Wh  '  ~2  72  Sln  2  Sin  2 


as 


-  f‘ 


(3.13) 


o 

This  solution  represents  a  fast  oscillation  with  frequency 
(  u>0+f)/2  being  modulated  by  a  slow  oscillation  (w0-f)/2. 
For  small  time  scales  t  <<  (wQ-f)-1,  and  for  a  frequency 
oo0  slightly  different  than  f,  an  approximation  to  (3.13) 
is 


w 


(us  +  f)t 

h~  (  4  t  Sin  - 2 - 


(3.14) 


an  oscillatory  solution  whose  component  grows  linearly  in 
time. 


3.1*2  Irrotational ,  Divergent  Wind  Field 


In  the  case  where  the  curl  of  the  wind  field  is 
zero,  only  the  divergence  terms  in  (3.6)  remain.  We  express 
the  wind  stress  in  the  form 


io  =  ( 0  ,<t>  0(y  )cosuj0t )  . 


(3.15) 


c  cos(w  t-<5)  +  <jsin(u)t-6) 

*hp  '  r,.  2  Zi  .2,2  .  .  5  2-il/2  ■  (3. IS) 


[<%  -  I*"'- 


!>4  *  (A2]1'2 


The  first  term,  proportional  to  cos(ui0t-  <5) ,  is  analogous  to 
(3.11a),  and  has  a  resonant  frequency  given  by  (3.12).  The 
second  term,  with  its  extra  factor  u>0,  has  a  different 
resonant  frequency 


OJ, 


(f2+c2)l/2  . 


(3.17) 
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there 


Fox  a  general  divergent,  rotational  wind  field, 
are  two  resonant  frequencies; 

u)r  =  (f2+c2)1/l2  ,  (3.18) 

where  the  plus  corresponds  to  the  first  two  inhomogeneous 
terms  in  (3.6),  while  the  minus  corresponds  to  the  third 
inhomogeneous  term.  So  we  have  two  resonant  frequencies, 
one  lower  and  one  higher  than  the  inertial  frequency.  This 
effect  may,  in  part,  be  responsible  for  the  two  near- 
inertial  peaks  resolved  by  Kundu  (1976)  in  his  kinetic 
energy  spectra  (see  Figure  2.3). 

3.2  INTERNAL  WAVE  RESONANCE  WITH  EKMAN  PUMPING 

In  Section  3.1  we  discussed  the  surface  layer, 
with  particular  emphasis  on  the  vertical  velocity  at  its 
base,  known  as  Ekman  pumping.  This  vertical  velocity  serves 
as  the  driving  boundary  condition  at  the  top  of  an  interior 
inviscid  layer.  We  assume  here  that  the  depth  of  the 
surface  layer  h  is  constant  in  space  and  time.  While  this 
assumption  is  not  generally  valid,  it  does  give  us  a  device 
for  decoupling  the  two  layers,  and  makes  the  problem  tract¬ 
able. 


We  start  with  a  set  of  interior  layer  equations, 
identical  with  (1.18)  except  that  we  set  the  stress  terms 
equal  to  zero.  The  interior  equations  can  be  combined  into 
a  single  equation  in  terms  of  vertical  velocity, 

9  t»zz  +  *2wzz  +  N2(z)(wxx+Wyy)  «  0  .  (3.19) 
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Ths,  boundary  conditions  are  w  «  0  at  the  bottom  z  =  -H , 
and  w  =  wjj  given  by  (3.8)  at  z  =  -h .  If  we  assume  a 
plane  wave  solution  proportional  to  exp  [  i(kx+  2,y+mz+  wt)  ] , 
then  we  get  a  relationship 

u>2  =  f2  +  N2(k2+  £2)/m2  .  (3.20) 

This  dispersion  relation  allows  a  free  internal  wave  to 
propagate  if  its  frequency  is  greater  than  the  inertial 
frequency  f.  This  means  that  w^c  in  (3.7)  does  not  con¬ 
tribute  to  the  production  of  free  internal  waves,  because 
whc  represents  oscillations  purely  at  the  inertial  frequency. 
On  the  other  hand,  the  second  term  for  w^p  in  (3.16)  does 
generate  free  internal  waves,  since  its  resonant  frequency 
is  greater  than  f. 

Those  frequency-wavenumber  spectral  components 
of  the  Ekman  pumping  that  satisfy  the  internal  wave  dis¬ 
persion  relation  are  said  to  be  in  resonance  with  internal 
waves.  We  can  visualize  these  resonances  in  terms  of  a 
modal  dispersion  relation 

o)j2(k,A)  =  f2  +  (k2+£2)N2H2/(jTr)2  ,  (3.21) 

where  we  have  replaced  m  in  (3.20)  with  j  tt/  H ,  and  have 
assumed  N  to  be  independent  of  depth.  Figure  3.1  is  a 
plot  of  (3.21)  for  the  first  four  modes  j  =  1,  2,  3,  4. 
To  be  concrete,  consider  a  simple  example  of  an  Ekman 
pumping  velocity  field,  whose  spatial  form  is  constant  as 
it  migrates  with  a  speed  U  in  the  x-direction,  and  is 
uniform  in  the  y-direction.  Then  we  could  write 

Wfc(x,t)  »  J  A(k)expik(x-Ut )dk  ,  (3.22) 
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Figure  3.1:  Dispersion  relation  for  lowest  four  modes 
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where  A(k)dk  is  the  amplitude  of  the  velocity  field  in 
a  wavenumber  band  dk  centered  about  a  wavenumber  k.  If  A(k) 
is  nonzero  at  a  wavenumber  which  satisfies  (3.15)  such  that 


o>j(k)  =  kU  ,  (3.23) 

then  an  internal  wave  is  excited  at  the  discrete  frequency 
oij.  Alternatively,  we  could  view  the  resonance  condition 
(3.23)  as  a  requirement  that  the  phase  speed  c  «  /k  of 
a  resonant  wave  must  be  equal  to  the  migration  speed  U.  A 
front  whose  length  scale  is  2Tr/k  ^  10$  m  must  have  a 

migration  speed  U  ^2.2  m/sec  to  excite  a  near-inertial  wave 
of  frequency  a,  1.4  x  10-4  sec-*. 

3.3  INTERNAL  WAVE  RESONANCE  WITH  THE  WIND  FIELD 


Miropol'sky  (1976)  investigated  the  generation 
of  internal  waves  using  a  two-layer  model.  He  derived  an 
expression  for  the  depth  dependent  horizontal  wavenumber 
spectrum  S(k,z,t)  of  vertical  velocities  in  terms  of  the 
spectral  tensor  of  wind-produced  surface  stress.  For  time 
scales  large  compared  to  the  integral  time  scale  of  the 
correlation  of  the  wind  stress,  the  expression  for  S(k,z,t) 
simplified  to  an  asymptotic  formula.  From  this  formula, 
Miropol'sky  concluded  that  the  resonant  growth  of  internal 
waves  occurs  only  when  the  spatial-temporal  spectra  of  the 
wind  stress  components  contain  a  set  of  frequencies  and 
wavenumbers  corresponding  to  the  dispersion  relations  of  the 
free  modes  of  the  internal  waves. 

Although  Miropol ’ sky' s  solutions  are  expressed  in 
analytic  form,  they  are  not  conducive  to  our  obtaining 
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physical  insight.  Instead,  we  review  some  results  from  a 
series  of  papers  by  Krauss  (1976a, b:  1978a, b;  1981)  which 
also  deal  with  the  resonant  production  of  near-inertial 
frequency  internal  waves  by  surface  stress.  In  these  papers 
he  presents  a  model,  a  method  of  numerical  solution,  and  a 
hierarchy  of  cases,  where  successive  cases  have  different 
assumptions  about  the  spatial  variability  of  the  wind 
stress.  Krauss  (1976a)  begins  with  a  set  of  equations  that, 
when  written  in  a  form  analogous  to  (1.18)  become 

2 

ut  +  f0v  =  -px  +  3  z(K3zu)  +  A7h  u  , 

2 

vt  -  fGu  =  -py  +  3  z(K3zv)  +  A  V  h  v  , 

wt  =  -Pz  +  *>  » 
b^  +  wN^  =  0  , 

ux  +  vy  +  wz  =  0  .  (3.24a-e) 

The  following  assumptions  are  made: 

1)  The  horizontal  eddy  viscosity  A  is  constant. 

2)  N  and  K  are  functions  of  depth  only. 

3)  The  viscous  terms  are  neglected  in  the  vertical  momentum 
and  buoyancy  equations. 


4)  Slip  is  allowed  at  the  sea  floor. 


We  note  that  the  hydrostatic  approximation  would  have  been 
valid,  but  our  discussion  in  connection  with  (1.16),  (1.17) 
showed  that  Krauss'  inclusion  of  the  wt  term  is  incon¬ 
sistent  . 

The  vertical  eddy  viscosity  K  =  (Ki;K2)  is  taken 
to  be  a  constant  in  an  upper  10  m  layer,  decreases 
linearly  from  Ki  to  K2  between  10  m  and  15  m,  and  remains 
constant  at  K2  from  15  m  down  to  the  bottom  (48  m  in  the 
Baltic  Sea).  The  equations  were  Fourier  transformed  in 
time  and  in  the  horizontal  direction,  and  were  solved 
numerically  in  depth. 

Krauss  expects  that  resonance  should  occur  at  the 
dispersion  lines  of  free  waves.  Figure  3.1  shows  a  repre¬ 
sentative  dispersion  relation  for  modes  1  through  4  in  an 
ocean  with  a  buoyancy  frequency  profile  that  is  independent 
of  depth.  The  dashed  line  represents  a  slice  at  a  single 
frequency  w  >  f.  Krauss  (1976b)  computed  equivalent 
vertical  displacements  along  such  a  slice,  at  different 
depths  and  for  different  eddy  viscosity  distributions 
( K 1 ; K2 ) -  Figure  3.2  shows  such  a  computation,  for  a 
divergent,  irrotational  wind  field.  Krauss  defined  the 
equivalent  vertical  displacement  to  be  the  amplitude  of  the 
vertical  displacement  if  the  particle  moved  as  w(z)eitt 
during  the  entire  period.  The  figure  caption  lists  addi¬ 
tional  details.  For  small  values  of  viscosity  (10;0.1), 
four  distinct  modes  are  produced.  As  the  viscosity  in¬ 
creases,  the  number  of  modes  decreases,  as  does  the  reso¬ 
nance  amplitudes.  When  the  viscosity  is  (100,100),  the 
internal  waves  are  entirely  suppressed. 


Before  proceeding  with  an  explanation  for  these 
results,  we  refer  the  reader  to  Figure  3.3,  which  shows  mean 
kinetic  energy  in  the  surface  layer  and  in  the  deep  layer, 
for  different  combinations  of  viscosity  (Ki;K2).  We  see 
that  for  (10;0.1)  and  ( 1 00 ; 1 ) ,  the  kinetic  energy  in  each  of 
the  layers  is  nearly  comparable.  As  the  ratio  K1/K2 
decreases,  the  resonances  broaden,  and  the  deep  layer 
kinetic  energy  decreases  relative  to  that  in  the  surface 
layer.  In  addition,  the  kinetic  energy  for  both  layers  in 
the  (100;1)  case  is  an  order  of  magnitude  weaker  than  in  the 
(10:0.1)  case. 

We  explain  these  results  as  follows.  A  certain 
amount  of  eddy  viscosity  is  necessary  to  transfer  momentum 
from  the  air  into  the  surface  layer.  If  the  Ekman  depth 
tt(  2Ki /f ) 1 /2  is  much  less  than  the  surface  layer  depth, 
then  there  is  no  stress  at  the  base  of  the  surface  layer. 
No  significant  damping  of  the  Ekman  pumping  occurs. 

The  smallest  value  Krauss  used  for  Kj  was  10  cm2/ 
sec  so  the  Ekman  depth,  12.6  m,  is  nearly  equal  to  the 
surface  layer  depth.  Significant  damping  does  occur.  From 
(3.16),  the  resonance  of  the  second  term  dominates  that 
of  the  first  term  for  small  values  of  c.  The  resonant 
amplitudes  of  the  first  and  second  terms  are  (2f)-1  and 
(2c)-1,  respectively.  As  c  increases,  the  resonant  am¬ 
plitude  of  the  second  term  decreases,  and  the  halfwidths 
of  both  terms  broaden.  This  effect  helps  to  explain  the 
resonance  broadening  in  Figures  3.2  as  3.3  as  Ki  increases 
from  10  to  100  cm2/sec. 

As  viscosity  increases  in  the  interior,  momentum 
is  allowed  to  diffuse  from  the  surface  layer  into  the 
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interior  at  a  faster  rate,  resulting  in  increased  damping 
of  the  Ekman  pumping.  From  Figure  3.2a,  we  see  that  the 
resonant  amplitudes  are  nearly  inversely  proportional  to  the 
mode  number.  Therefore  as  the  interior  viscosity  K2 
increases,  the  reduction  in  resonant  amplitudes  in  the 
interior  is  nearly  independent  of  mode  number. 

Krauss'  (1978a)  interpretation  of  these  figures 
is  of  value.  We  can  think  of  these  diagrams  as  being 
response  functions  Hq(k,  w)  such  that 

q(k,u>)  =  Hq(k,  ui)T0(k,u))  , 

for  any  physical  variable  q.  The  resonance  peaks  correspond 
to  free  internal  waves.  Typical  gain  factors  range  from  2 
to  10.  Surface  stress  components  whose  wavenumbers  fall 
outside  the  resonant  peaks  generate  forced,  or  trapped 
internal  waves,  in  the  sense  that  they  are  not  free  to 
propagate  away  from  their  local  source.  Stress  components 
whose  wavenumbers  fall  in  the  vicinity  of  resonant  peaks 
generate  freely  propagating  internal  waves.  This  is  the 
subject  of  the  next  section. 

3.4  PROPAGATION  OF  NEAR- INERTIAL  FREQUENCY 

INTERNAL  WAVES 

3.4.1  WKB  Solution 

We  wish  to  analyze  the  propagation  of  low  fre¬ 
quency  internal  waves.  We  start  with  an  Ekman  pumping 
vertical  velocity  at  the  base  of  the  surface  layer  as 
the  source  of  internal  waves,  and  examine  their  propagation 
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downward  into  the  interior.  We  assume  that  the  vertical 
velocity  at  z  =  -h  is  of  the  form: 

w(z  =  -b)  =  Wjj(x,t)  exp[i(kx  -  ut)],  (3.25) 

where  the  amplitude  Wh(x,t)  is  a  weak  function  of  x  and  t; 
3xWh/Wji  and  3twh/Wh  are  much  less  than  k  and  w  t  respec¬ 
tively  . 


The  equation  of  motion  from  (3.19)  is 

9z(3t  +  f2)W  +  n2(z>(3x  +  9y)w  =  0  (3.26) 

With  no  loss  of  generality,  we  let  the  horizontal  component 
of  propagation  be  in  the  x-direction.  Then  we  can  assume  a 
solution  for  w  of  the  form 

w(x,z,t)  =  wo(x,z,t)exp£iS(x,z,t)[]  .  (3.27) 

S(x,z,t)  is  a  phase  function  which  varies  much  more  rapidly 
with  x,  z,  and  t  than  the  amplitude  function  wQ.  We  let 


3  XS  *  k  , 

3  ZS  =  m  , 

3tS  a  -  0)  ,  (3.28a-c) 

so  that  (3.16)  becomes 

32zwD  +  m2w0  »  0, 

m2(z)  a  k2N2( z)/(w2  -  f2)  .  (3.29a,b) 


3-18 


In  order  that  m  be  real  valued,  we  require  w  >  f  .  The 
terms  k  and  to  are  the  horizontal  wavenumber  and  frequency, 
and  m(z)  acts  as  a  local  vertical  wavenumber. 

If  N(z)  is  a  smoothly  varying  function,  that  is 
N'(z)/N(z)  <<  m(z),  then  the  WKB  approximate  solution  is 
appropriate; 

w  =  A  m(z)”1/2  exp£i(kx  + J  mdz  -  wt)J  .  (3.30) 

Here,  the  integration  limits  are  -h  and  z,  and  A  is  a  con¬ 
stant  to  be  determined  by  the  boundary  condition  at  z  *  -h. 
For  concreteness,  we  set 

A  -  Wh  raj,1/2,  (3.31) 


where  m^  *  m(z  *  -h),  so  that  (3.30)  agrees  with  the  bound¬ 
ary  condition — and  since  we  are  concerned  with  only  downward 
propagation  we  retain  only  the  lower  minus  sign  in  (3.30). 
Then  the  real  valued  solutions  for  u,  v,  and  w  (see  Kundu; 


1976)  are 


cos4> 


=  W  - 
h  w 


(mhm) 


sin<|> 


(3.32a-d) 
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3.4.2 


Energy  Piupagation 


Energy  propagates  downward,  away  from  points  on 
the  plane  z  *  -h ,  along  characteristic  curves.  The  propaga¬ 
tion  proceeds  at  the  local  group  velocity  cg,  given  by  Cg  = 
(3w/3k,  3w/3m).  Using  the  characteristic  equation  (3.19b), 
we  compute  the  group  velocity  to  be 


2  .2 

cg  =  (1,  -(u>2-f2 )1/2/N( z)  )  .  (3.33) 

For  near  inertial  frequencies,  tu-f  <<  f  holds,  and  the 
magnitude  of  the  group  velocity  is  approximately 

l-g 1  ~  *\w-f)/k  ,  (3.34) 

and  if  N(z)  >>  f,  then  the  angle  0  of  the  group  velocity 
vector  with  respect  to  the  horizontal  is  very  small: 


0  —  tan 


-l/y^m-f  )f  K  /2  ( u)-f  )f 


N(z) 


N(  z) 


(3.35) 


Pollard's  (1980)  analysis  of  current  meter  records 
give  us  some  insight  into  the  propagation  of  near-inertial 
internal  waves.  Figure  3.4  illustrates  horizontal  and 
vertical  coherences,  and  vertical  phase  differences  among 
three  moorings,  separated  horizontally  by  about  50  km. 
Figure  3.4c  shows  statistically  that  significant  values  of 
vertical  coherence  correspond  to  upward  phase  propagation; 
phase  differences  all  lie  between  20°  and  70°.  The  vertical 
component  of  group  velocity  is  therefore  directed  downward, 
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Figure  3.4:  Coherence  and  phase  of  frequencies  between 
(20.89)-!  and  (16.25)-!  cph.  For  each  pair 
of  vector  time  series,  four  values  of  coherence 
and  phase  can  be  calculated  (Webster,  1968). 
Each  vertical  bar  joins  the  minimum  and 
maximum  of  the  four  values.  Horizontal  coher¬ 
ence  values  (a)  are  given  for  pairs  of  records 
at  the  same  depth  on  different  moorings. 
Vertical  coherence  (b)  and  phase  (c)  are  given 
for  pairs  of  records  on  the  same  mooring  20m 
apart.  Phases  are  modified  by  +  90°  (cf. 
Webster,  1968)  to  give  the  angle  by  which  the 
second  series  leads  the  first  series.  The  95% 
limit  for  the  coherence  is  0.45.  From  Pollard 
(1980) . 


auu  the  inertial  oscillations  at  32,  52,  and  72m  are  wind- 
forced  oscillations  dispersing  away  from  the  surface.  From 
measurements  of  horizontal  and  vertical  phase  differences, 
Pollard  estimated  the  group  velocity  to  be  on  the  order  of 

Cg  '■o  (1-20  km/day,  0.03-3  m/day), 

and  suggests  these  small  values  are  responsible  for  the 
rapid  loss  of  horizontal  coherence  with  increasing  depth  in 
Figure  3.4a.  Using  the  maximum  group  velocity  estimates,  it 
takes  6-7  days  for  near-inertial  energy  to  propagate  20m 
vertically,  and  3-4  days  to  travel  from  one  mooring  to 
another.  During  these  several-day  periods  and  over  tens  of 
kilometers  the  variability  in  the  upper  ocean  reduces  the 
coherence  of  the  inertial  motions.  The  group  velocity  is 
sensitive  to  horizontal  variations  in  stratification,  and 
coherence  may  be  destroyed  by  a  number  of  processes,  includ¬ 
ing  Kelvin  Helmholtz  instability,  turbulence,  and  mean 
horizontal  as  well  as  vertical  shear. 

3.4.3  Vertical  Shear 


Equations  (3.32b,c)  describe  a  horizontal  velocity 
vector  that,  at  any  given  point  in  space,  rotates  in  a 
clockwise  ellipse  in  time.  The  major  axis  is  in  the  direc¬ 
tion  of  propagation,  and  the  minor  axis  is  smaller  by  the 
ratio  f/w. 


The  horizontal  velocity  vector  turns  as  a  function 
of  depth.  Looking  down,  the  turning  is  clockwise  with 
increasing  depth  for  upward  propagation,  and  counterclock¬ 
wise  for  downward  propagation.  Kundu  (1976)  computed  the 
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turning  rate,  but  we  suspect  there  is  a  misprint  in  his 
text.  The  corrected  computation  follows  below.  Define  <j>  as 
the  angle  between  the  horizontal  velocity  vector  and  the 
direction  of  horizontal  propagation,  the  x-axis.  Then  from 
(3.32b,c) 

tan<j>  -  ~  —  tan(kx  -  /  mdz  -  ust ) 

u  w  /  (3.36) 


The  turning  rate  can  be  written 


dd> 

dz 


1  +  tan^(kx  -  wt  -  J mdz) 


1  +  tan 
_  1  +  (f/w) 


tan2(kx  -  cat  -  J  mdz 


) 


(3.37) 


With  the  approximation  u  s  f,  (3.37)  becomes 

:=  -m(z)  ,  (3.38) 

which  simply  states  that  the  velocity  vector  turns  by  about 
360°  in  a  distance  equal  to  the  vertical  wavelength. 

We  can  also  compute  the  vertical  shear; 

s2  =  uz  +  vz  -  (m3  +  mz/4m>  *  (3.39) 

k 

Prom  the  WKB  approximation  Nz/N  <<  m,  we  are  justified  in 
neglecting  the  second  term  in  parentheses  in  (3.39).  Then 
combining  (3.29b)  and  (3.39)  gives 
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(3.40) 


(z) - 


Wh  rahkN3^2> 


Here  we  have  derived  a  standard  WKB  result,  that  S2  is 
proportional  to  We  can  interpret  these  results  in  a 
simple  manner.  From  (3.22b),  we  could  express  the  amplitude 
of  the  horizontal  current  as 


|U|  =  Vfh(mhm)1/2  k"1 


(3.41) 


and  the  shear  becomes 

S(z)  =  j  U(z) 1  m(z) 


(3.42) 


In  the  case  of  a  uniformly  stratified  ocean,  the  vertical 
wavelength  Xv  is  2iT/m,  and  (3.32)  might  be  written 

S  =  2tt|U|  /Av  .  (3.43) 

We  emphasize  that  these  results  apply  for  a  monochromatic 
internal  wave  that  would  accompany  an  ideal  sinusoidal  form 
of  Ekman  suction  (3.15),  hence  an  ideal  sinusoidal  surface 
stress. 
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Section  4 

SUMMARY  AND  FUTURE  RESEARCH 


Most  of  the  research  to  date  deals  with  the 
resonant  excitation  of  internal  waves  by  wind  forcing.  The 
emphasis  has  been  on  the  growth  rate  of  resonant  vertical 
modes,  where  the  internal  wave  field  responds  to  a  steady 
state  wind  field.  A  steady  state  wind  field,  though,  is  not 
a  realistic  model  of  the  atmosphere.  A  wind  forcing  event 
with  a  horizontal  wavenumber  component  k  and  propagation 
speed  U  may  be  in  resonance  with  the  internal  wave  field 
for  a  short  while.  However,  if  U  is  not  constant,  the 
response  will  quickly  slip  out  of  that  resonance,  and 
perhaps  into  a  different  one. 

For  these  reasons,  it  appears  that  near-inertial 
frequency  internal  waves  are  transient  responses  to  variable 
wind  forcing.  An  important  problem  that  should  be  addressed 
is  the  response  due  to  wind  events  that  propagate  with 
variable  speed,  or  whose  structure  or  length  scales  evolve 
in  time.  A  particular  aspect  of  this  problem  would  be  to 
determine  the  degree  of  statistical  stationarity  required  of 
a  wind  field,  in  order  to  elicit  a  resonant  response. 

To  be  able  to  study  the  transient  response  of  the 
internal  wave  field,  we  are  developing  a  numerical  model  of 
a  wind-forced,  viscous,  stratified  ocean.  The  model  equa¬ 
tions  are  similar  to  those  of  Krauss,  but  the  method  of 
solution  is  different.  By  employing  a  finite  difference 
scheme  which  marches  the  solutions  forward  in  time,  it  is 
possible  to  study  the  generation  of  near-inertial  motions 
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and  shear.  The  model  does  not  impose  severe  'imitations  on 
the  functional  form  of  the  wind  field,  and  does  not  require 
that  it  be  periodic  in  time.  Details  concerning  the  model 
and  results  will  be  given  in  Rubenstein  (1981). 
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APPENDIX  A 


We  solve  the  equation 

Ut  +  if  U  -  KUZZ  =  Q(z.t)  +  U0(z)6(t)  ( A .  1 ) 

for  a  complex  source  distribution  of  momentum  Q(z,t) 
and  an  initial  complex  velocity  LJ0  =  U(z,0).  With  a 
horizontally  homogeneous  source  distribution,  we  let 
Q  =  2i0(t) 6(z),  where  the  factor  of  two  enters  because  we 
have  extended  our  domain,  in  a  virtual  sense,  above  the 
plane  z  =  0.  We  then  write 

Ut  ♦  If  0  -  KVZZ  -  2To(t)6(Z)  .  ( A .2) 

We  apply  a  Laplace  transform  operation  to  (A. 2)  with  respect 
to  t :  j 

sU(z,s)  -  Uq( z )  +  ifU(z.s)  -  KUzz(z,s)  =  2tq(s)6(z)  (a. 3) 

Next  we  take  the  inverse  Fourier  transform  of  (A. 3)  with 
respect  to  z; 

sU(B,s)  +  (if  +  4t2S2  K)U(6  ,s)  =  U  (6)  +  2t0(s> 

We  can  solve  for  U(g,s) 

U(  8 ,  s )  =  [s  +  (if  +  4n282K)]  1  [uq(3)  +  2tq(s)] 

and  take  the  inverse  Laplace  transform  to  get 


(A. 4) 


•  (A. 5) 


U(  6  ,  t )  =  exp£-(if  +  4772S2K)tJ  UQ(S) 

+  /  exp  T-  (  i  f  +  47T282K)tll  *  2X  (t) 

'  '  (t)  ° 


(A. 6) 


Here  the  asterisk  indicates  the  operation  of  convolution 
with  respect  to  t.  Now  we  take  the  Fourier  transform  of 
(A. 6)  with  respect  to  z;  we  use  the  fact  that  the  Fourier 
transform  of  exp  (-Tra  g2)  is  a-l/2  exp(-Trz2/a)  and  the 
convolution  theorem  to  yield 

U(z,t)  =  e"lft(4'aKt)~1'/2  [exp(-z2/4Kt )]  *  U  (z) 

(z) 

+  re_ift(4-fiKt  )_1//2  exp(-z2/4Kt )"]  *  2t  (t)  (4  7) 

J  ( t )  0 

In  the  special  case  of  a  steady  state  situation, 
we  set  Ut  =  0  in  (A. 2)  and  Fourier  transform  with  respect 
to  z; 


U(6)  = 


2  2 

if  +  4tt 


(A. 8) 


We  let  a2  =  if/K,  so  that  (A. 8)  may  be  written 


U(B)  =  —  2a 

Ka  (  2tt 6  )2  ‘ 


The  inverse  Fourier  transform  of  (A. 9)  is 
—  i  7i  /  4 

T  6  ' 

U(z)  =  - — -  exp  f-  ( 1  +  i)(f/2K)1/2  | 

(fK)1/2  L 


(A. 9) 


(A. 10) 
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This  result  is  the  well  known  Ekman  spiral.  At  the  surface, 
with  z  =  0,  the  phase  e-*'*/4  indicates  that  the  current  is 
45°  to  the  right  of  the  wind  stress. 


Parenthetically,  it  is  helpful  here  to  define  the 
convolutions  in  (A. 7).  With  the  appropriate  symmetry  of 
u0(z),  the  convolution  with  respect  to  z  is  defined  by 


/ 


f ( z )  *  UQ( z)  =  /  f (y)u0(z-y)dy 


(A. 11) 


Also,  the  convolution  with  respect  to  t  is  defined  by 

t 


f(t)  *  T 


o(t)  -  f 

•Jn 


f  (y)'ro(t~y)dy 


(A. 12) 
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